Unified Explanation of the Solar and Atmospheric neutrino Puzzles in a 

supersymmetric SO(10) model 
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It was recently suggested that in a class of supersymmetric SO(10) models with Higgs multiplets 
in 10, and a single 126 + 126 representations, if the 126 contributes both to the right handed 
neutrino masses as well as to the charged fermion masses, one can have a complete prediction of the 
neutrino masses and mixings. It turns out that if one chooses only one 10, there are no regions in the 
parameter space where one can have a large — v T mixing angle necessary to solve the atmospheric 
neutrino deficit while at the same time solving the solar neutrino puzzle via the z/ e <-> oscillation. 
We show that this problem can be solved in a particular class of SO(10) models with a pair of 10 
multiplets if we include the additional left-handed triplet contribution to the light neutrino mass 
matrix. This model cannot reproduce the mass and mixing parameters required to explain the 
LSND observations neither does it have have a neutrino hot dark matter. 
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Strong indications in favor of non-vanishing neutrino 
masses are emerging from several experiments: (i) the 
deficits of solar neutrino flux observed by the four so- 
lar neutrino experiments Homestake, Kamiokande, Su- 
perKamiokande, SAGE and GALLEX fy] compared to 
the standard solar model calculations |j can be under- 
stood if neutrinos are massive and the electron neutrinos 
emitted by the sun oscillate to another neutrino species, 
and (ii) the atmospheric muon neutrino deficits observed 
earlier by Kamiokande, 1MB and Soudan II H experi- 
ments and confirmed recently by Super Kamiokande can 
be understood if oscillates similarly. The LSND ||] 
results have provided the first laboratory indication of 
v & oscillation and if confirmed by KARMEN ||, 
would seal the case for non-zero neutrino masses, in an 
unequivocal manner. 

As is well known, the solar neutrino deficit can 
be explained in terms of the matter induced resonant 
Mikheyev-Smirnov-Wolfenstein (MSW) oscillation |j) for 
two choices of masses and mixing angles M. Our in- 
terest here is in the so called small angle solution for 
which A? 
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1.0) x 10~ 5 eV 2 and 2 x 10~ 3 < 
e(1 n i A !„ , The atmospheric neutrino deficit 
could be due to either v T or v e oscillation. 

Preliminary indications from the electron energy distri- 
bution in SuperKamiokande favors <-> v T oscillation. 
Similarly a preliminary fit to all the atmospheric neutrino 
data (sub-GeV, multi-GeV including the zenith angle de- 
pendence) seems to require 2 x 10 -4 < m 2 T (eV 2 ) < 10~ 2 



with sin 29 b 



0.6-1.0 



Note the hierarchical pat- 



tern of mass differences. The LSND results require that 
0.3 eV 2 < Am 2 M < 10 eV 2 with the mixing angle in the 
few percent range. If we accept the above results, it is 
clear that with only three neutrinos, it is not possible 
to explain the three results (i.e. solar, atmospheric and 
LSND) simultaneously Therefore within conventional 
grand unified theories with three generations, one may 
hope to understand only two of the above results. Fur- 
thermore, since hierarchical mass patterns for neutrinos 



is a generic feature of theories that implement the see-saw 
mechanism Q, the solar and the atmospheric neutrino 
data appear more amenable to theoretical understanding 
in simple models. 

It is the goal of this paper to present a simple grand 
unified scheme (GUT) that leads to the supersymmetric 
standard model at low energies and predict Am 2 and 
sin 2 28 values in the above range for v e <-> and «-> v T 
sectors so that we have a theoretical understanding of the 
solar and the atmospheric neutrino data. We believe this 
result to be significant since we do not use any extra 
fcrmions nor any extra symmetries for the purpose. 

The simplest GUT theory that leads naturally to small 
neutrino masses via the see-saw mechanism || is the 
SO(10) model where the local B-L symmetry is broken 
by the 126 + 126 representation. It also has another at- 
tractive feature that it leads to automatic R-parity con- 
servation so that unwanted (and uncontrolled) baryon vi- 
olating interactions of the MSSM are forbidden and one 
obtains a stable LSP which can act as the cold dark mat- 
ter of the universe. The minimal set of Higgs multiplets 
needed to break all gauge symmetries of the theory while 
keeping supersymmetry unbroken down to the weak scale 
is: 45 + 54 (denoted by A and S) 126 + 126 (denoted 
by A and A) and a single 10, denoted by H. 

It was shown, sometime ago flo| , |ll| that in this mini- 
mal model, all Yukawa couplings and Higgs vevs respon- 
sible for fermion masses and mixings (a total of twelve 
parameters in all in the absence of CP-violation) are com- 
pletely determined by the quark and lepton masses and 
the quark CKM angles. As a result the light and heavy 
Majorana mass matrices for the neutrinos are completely 
determined except for the overall scale Vr, the scale of 
B-L symmetry breaking, provided one assumes the sim- 
ple see-saw formula (to be called type-I see-saw formula 
§) 
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This enables a complete prediction of neutrino mixing 
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angles and any two neutrino mass ratios. Every choice 
for the signs of the various charged fermion masses lead to 
distinct scenarios and separate predictions. It was found 
that there were predictions that could accommodate only 
small angle MSW solution to the solar neutrino puzzle 
but not the atmospheric neutrino puzzle. The reason 
was that the maximum value for sin 2 20 MT mixing angle 
predicted by this model was less than 0.3 or so, where 
as the present 99% confidence level fits seem to require 
sin 2 20 MT ~ 0.60 or higher g. 

One may try to take advantage of the fact that in most 
left-right and SO(10) models, a generalized see-saw for- 
mula for neutrino masses holds |L2] (to be called the type 
II see-saw formula): 

M v = fvi - M° M-i \M° ] T , (2) 

v 2 

(where vl — an d is induced as long as there are 

54 dimensional Higgs multiplets in the theory) and see if 
it is possible to obtain larger values for sin 2 29 llT . Such 
models would have two free parameters, vr and vl- In 
SUSY GUT models, coupling constant unification includ- 
ing threshold corrections puts vr from 10 13 to 10 15 GeV 
range leading to vl — 10~ 2 to 1 eV for A ~ 1. The 
two terms then give comparable contributions and we 
have two parameters {vl,vr} that determine the neu- 
trino masses and mixings. We have made an extensive 
numerical analysis of the predictions of this model for 
neutrino masses and were unable to find any reasonable 
values of vl and vr which can accommodate both the 
small angle MSW solution to the solar neutrino problem 
as well as the <-> v r oscillation solution to the atmo- 
spheric neutrino puzzle. Thus neutrino experiments may 
play the role similar to the role that proton decay exper- 
iments played in ruling out minimal non-supersymmetric 
SU(5) model. 

We therefore are led to consider a slight generalization 
of the above idea and consider an SO(10) model with two 
10-dim. Higgs multiplets instead of one as in the minimal 
model. The rest of the field content is the same. The 
remainder of the paper will be devoted to studying the 
neutrino masses in this model. The low energy theory 
in this model is the MSSM with the Higgs doublets in 
general being linear combinations of the doublets in the 
10's (denoted -Hi, 2) and the 126. We will assume the 
following specific form for them. 

H u = ai H u {lQx) + a 2 H u {126) + a 3 H u {126) 

H d = (3 x H d {W 2 ) + /3 2 #d(T26) + f3 z H d {12Q) (3) 

How the light doublets arise with this specific form is of 
course related to the difficult problem of doublet-triplet 
splitting in SO(10) models which is not addressed here. 

Let us now discuss how the neutrino mixing angles can 
be extracted from this model. The first point is that the 
most general Yukawa superpotential of the model given 
by, 



Wy = h it ab tpa Hi + f ab ip a ip b A, (4) 

where ip a (a — 1, 3) represent the 16 dimensional spinors 
corresponding to the three family of fermions. Since 
SO(10) symmetry implies that hi and / are symmetric 
matrices, (we ignore CP violation from Yukawa sector), 
we can diagonalize any one of them and we have fifteen 
free Yukawa coupling parameters in terms of which the 
fermion masses and mixings are expressed as follows: 

M u = hi v u + f k u M d = h 2 v d + f K d 

M v d =hv u - 3/ k u Mi = h 2 v d - 3/ n d (5) 

Using these relations, we find that at the GUT scale, we 
have 

M v d = n{Mi - M d ) + M u (6) 
M„ LVL = r 2 (M d - Mi) (7) 
M Nr = r 3 (M d - Mi) (8) 

where n = ^;r 2 = ^ and r 3 = ^ and M^ is 
assumed to denote the fi)L contribution to the neutrino 
mass matrix. From the above equations it is clear that we 
need to supply six parameters to determine the neutrino 
masses and mixings and they are the three miximg angles 
in the charged lepton mass matrix and (i — 1,2,3). We 
demand that the three charged lepton mixing angles are 
zero. We then scan the parameter space for r.; to see if 
any desirable solution exists. 

To proceed with this program, first note that the above 
relations between fermion masses hold at the GUT scale. 
So, we extrapolate the observed values of quark and lep- 
ton masses to the GUT scale, using simple analytic for- 
mulae given by Naculich |Q . We work in a basis where 
M u is diagonal and M d = V ckm D d V^ km . At the GUT 
scale the diagonalized values for the masses in GeV and 
the values of the angles are 

m u = 0.0011 m c = -0.3785 m t = -112.34 

m d = 0.00131 m s = 0.0148 m b = -1.177 

m e = 0.0003 to m = -0.0699 m T = 1.183 

8x2 = -0.2210 S13 = 0.0040 s 23 = 0.03 1 

where s\ 2 is the Cabbibo angle, S13 and S23 are roughly 

the V u b and V cb elements of V c km- 

In the basis we are working, Mi is diagonal. Further- 
more, since the signs of the fermion masses are arbitrary, 
we choose a basis where the various fermion masses have 
the signs as given above. We then use Eq.(|^) and Eq.(j|) 
for each of the cases, to obtain the neutrino masses and 
mixing angles. 

Note that we still need to know vl = ^r 2 n d and 
Vr = Ar 3 n d . One can use theoretical arguments for the 
orders of magnitude of the parameters vl and vr, that are 
plausible. Note, for instance that since the value of the 
induced vev vl — f^, for vr in the range of 10 13 — 10 16 
GeV, vl — 1 — 10~ 2 eV is quite reasonable. One way 
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to determine vr is to use the unification constraint as it 
applies to the minimal model. We assume that the the- 
ory below the Vr scale is the minimal supersymmetric 
standard model (MSSM). Since different choices of the 
particle spectrum above the intermediate scale give dif- 
ferent values of vr, we use another method to constrain 
this parameter. 

Baryogenesis Constraints on the Scale vr 

A very simple mechanism for baryogenesis in SO(10) 
models is to generate a lepton asymmetry at a high tem- 
perature via the decays of the right-handed Majorana 
neutrinos and have this lepton asymmetry converted to 
a baryon asymmetry [jl5j by the sphaleron processes. An 
important necessary condition for this to happen is that 
at-least one of the right handed Majorana neutrinos must 
have a decay rate that is slower than the expansion rate 
of the universe when T ~ Mj\r. The general formulae are: 
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r 



J2b ^i,ab 



8n 



M 2 

M Na < 1.73 (g*) 1/2 -rf±- (9) 
Mpi 



Since in our model the Yukawa couplings arc all predicted 
in terms of vev's v u ,Vd and Kd, we can obtain a lower 
bound on Mff t if we know the vev's and using the pre- 
dicted value for / matrix, we can then deduce vr. Since 
in our analysis is independent of v u and Vd, the only con- 
straint on them is that \Jv u + v\ + k u + k 2 , = 246 GeV. 
Using the fact that we have chosen 7*1 = 40.3 (see be- 
low), we can get Kd — 1 GeV. Using them, we find that 
v R > 10 14 GeV. 

Prediction for Neutrino masses and mixings 

Using Mathematica we have scanned over all possible 
choices for the signs of the charged fermions to see if 
there is a prediction that fits the requirements of both 
the Solar and the atmospheric neutrino puzzles. In Fig- 
ures (1) and (2) we we plot the A m 2 ea and sin 2 2Q eu as 
functions of Vl and vr. In Figure (3) and (4) we plot 
A tti 2 r and sin 2 20 MT as functions of vl and vr. Cases 
{A,B,C} have r 3 = {1.81, 1.89, 1.98} x 10 13 respectively 
We see that sin 2 28 Mr is the most sensitive function of 
vl and we find acceptable solutions, displayed in Figure 
(4): {m Ve , m Ufi , m VT \ 
eV and 



-{0.063,3.087,10.88} 10~ 3 in 



-0.989 -0.081 -0.123' 
i T „ = | -0.147 0.539 0.829 
0.001 -0.838 0.545 



(10) 



Which gives sin 2 2Q ea = 2.8 x 10~ 2 and sin 2 26,, T = 0. 
and the Yukawa couplings hi/ 2 and / are, 



hi = — 

V u 



-0.064 0.123 0.181 
0.123 -3.721 1.488 
0.181 1.488 -17.285 



(11) 



1000 k d 




-0.761 
20.735 
-9.230 



-1.126 
-9.230 
-589.692 



(13) 
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Where we had r\ = 40.3, r 2 = 3.15 x 10"", r 3 
1.89 x 10 13 . Using this, explicit determination of hi/ 2 
in combination with baryogenesis constraint Eqn. (^) , we 
obtain the lower limit on vr > 10 14 GeV as stated ear- 
lier. A few comments are in order on other aspects of the 
SUSY SO(10) model characterized by the superpotential 
in Eq.(|) 

(i) The doublet-triplet splitting in this model has the 
non-trivial property that it leads to realistic fermion mass 
spectrum in contrast with the Dimopoulos-Wilczck (DW) 
mechanism. The point is that in the DW case the MSSM 
doublets arise from 10 dimensional SO (10) multiplets 
thereby leading to incorrect mass relations — = — , 

J Trip Tn s 5 

which is off by a factor 10 or so. In contrast, in our 
model the low energy MSSM doublets are admixtures of 
doublets in 10 and 126 and is therefore free of such dif- 
ficulties. 

(ii) It is also worth emphasizing that, the near maximal 



mixing angle for 



sectors needed to explain the 



atmospheric neutrino data is very hard to obtain with 
the type I see-saw formula as has been clear in many 
studies [ 10 |ll|,|lj| . One generally needs heavier vector like 
quarks [16 1 for this purpose. Thus our analysis would 
speak in favour of the type II see-saw formula which puts 
constraints on the SO(10) model building, 
(iii) Strictly speaking the prediction of neutrino masses 
is sensitive to the renormalization of the see-saw formula 
Jl7| . However, in our case, the Yukawa couplings are 
so small that (as can be seen from Eqs. 11, 12 and 13 
) this extrapolation does not noticeably alter the above 
predictions at low energies. 

In conclusion, we have shown that the use of a type 
II see-saw formula in a next to minimal SUSY SO(10) 
model without extra matter multiplets or extra symme- 
tries can explain both the solar and atmospheric neutrino 
deficits but not the LSND results. Thus if the LSND re- 
sults are confirmed by KARMEN experiment, this class 
of SUSY SO(10) models (minimal and next to minimal) 
cannot accomodate it simultaneously with the solar and 
the atmospheric neutrino results and alternative theoret- 
ical frameworks must be investigated. 
This work was supported by the National Science Foun- 
dation grant number. PHY-9421386. We are grateful to 
Bhaskar Dutta for critical comments on an earlier version 
of the paper. 
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„ for various Vl and vr. 



Predicted sin 2 20 e/J for various vl and vr. 
Predicted AmJ T for various vl and Vr. 
Predicted sin 2 20 Mr for various vl and vr. 



